Efimov states embedded in the three-body continuum 
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We consider a multichannel generalization of the Fermi pseudopotential to model low-energy 
atom-atom interactions near a magnetically tunable Feshbach resonance, and calculate the adiabatic 
hyperspherical potential curves for a system of three such interacting atoms. In particular, our 
model suggests the existence of a series of quasi-bound Efimov states attached to excited three- 
body thresholds, far above open channel collision energies. We discuss the conditions under which 
such states may be supported, and identify which interaction parameters limit the lifetime of these 
states. We speculate that it may be possible to observe these states using spectroscopic methods, 
perhaps allowing for the measurement of multiple Efimov resonances for the first time. 



Three resonantly interacting particles may form long- 
range bound states, called Efimov states even if 
the short-range interparticle interaction supports no two- 
body bound states. After years of tantalizing work on 
atomic "^He trimers and halo nuclei (see and references 
therein) , the first strong experimental evidence of Efimov 
physics was recently found in an ultracold thermal gas of 
cesium atoms 0]. This observation was made possible 
by utilizing a magnetically tunable Feshbach resonance 
to precisely control the two-body scattering length over a 
substantial range, and observing a resonant feature in the 
atom-loss rate due to three-body recombination aaa. 
Since hyperspherical studies have provided important in- 
sights in the past @, 0, Sl j '^^ ^i'^ motivated in this Let- 
ter to investigate the nature of universality within the 
adiabatic hyperspherical representation for multichannel 
two-body interactions, obtaining a complete view of the 
complicated energy landscape. 

Multichannel systems have a considerably richer struc- 
ture than their single-channel counterparts admit- 
ting multiple three-body thresholds, quasi-bound two- 
body channels, and multiple length scales. In this Let- 
ter, we present a Green's function method for obtaining 
the adiabatic potential curves, and discuss some immedi- 
ate consequences. Note that the zero-range multichannel 
method developed initially by Macek and Kartavtsev to 
study three-body recombination near a narrow two-body 
Feshbach resonance is the closest existing study to 
the present analysis, but it did not consider our main 
point of focus, namely the Efimov effect for excited three- 
body thresholds. 

Efimov states attached to excited three-body thresh- 
olds are absent in single channel models and difficult to 
infer from other multichannel approaches |12| . For these 
excited thresholds, at least one of the three atoms re- 
sides in an excited Zeeman or hyperfine state. We label 
these one-atom states |1) and |2) with energies -Ei = 
and E2 = e, respectively. In particular, we show that 
under the rather general stipulation that a quasibound 
two-body state is degenerate with an excited two-body 
threshold (see Fig. [1]), and that the coupling between 
product states |11) and |22) is much smaller than all other 
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FIG. 1; (color online) A schematic representation of our two- 
body model is shown. As the range ro is taken to zero, the 
potential is regularized by the derivative in Eq. ([1} 



couplings, a series of Efimov states should emerge that is 
attached to an excited three-body threshold. These Efi- 
mov states are of a unique character and can be viewed 
as fully three-body Fano-Feshbach resonances embedded 
within the three-body continuum far above the energy 
of open channel collisions, in contrast with the single- 
channel Efimov resonances predicted and recently ob- 
served that were interpreted as shape resonances [3, 01 ■ 
Since the interatomic scattering length for ground-state 
atoms will not in general be resonant, the gas will be 
comparatively stable with respect to the scaling law 
for three-body recombination @, 0, @l • Efimov states at- 
tached to the excited threshold could be directly accessed 
via photoassociation and observed through the measure- 
ment of atom-loss rates, thus opening a new toolkit of 
powerful spectroscopic techniques [14 1. 

As illustrated in Fig. [l] the two-body system has a 
set of product states {\(t = 1),|(t = 2),|(t = 3)} = 
{|11),(|12) + |21))/V2, |22)} with respective energies 
{Ea-} — {0,e, 2e}. The zero-range two-body potential 
is now written in the {\cr)} basis as. 
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This potential is the natural extension of the regular- 
ized Fermi pseudopotential which imposes the Bethe- 
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Peierls boundary condition on the two-body wavefunc- 
tion, ■i/;(r) C(l — a/r) as r ^ 0. When all channels 
are energetically open, the matrix A is simply minus the 
scaled two-body reaction matrix ifg. The physical reac- 
tion matrix is obtained from the multichannel quan tum 
defect theory (MQDT) channel closing formula jlSl. 
Note that unlike other multichannel treatments that re- 
tain only one open background channel and a one reso- 
nant molecule channel I4I, we label our channels explic- 
itly by the internal states of the atoms, and our Hamilto- 
nian allows for scattering at excited thresholds between 
atoms with different internal states. 

Unless explicitly stated, we choose the energy inde- 
pendent symmetric matrix to have elements Aij such 
that Ai3 0, and A33 An ^ 3.75 ao giving e/h — 
h"^ / (hmA^^) ~ 3 GHz and placing a quasi-bound state at 
E = e. The quantity ^33 is the length scale by which we 
measure all other quantities for the remainder of this Let- 
ter. We choose A12 = A23 = A33/4 to give a reasonable 
resonance width Ab ^ — 8G, and A22 = 2A33 in order 
to illustrate the effect of a quasi-bound two-body state 
on the three-body potential energy curves. These num- 
bers are chosen to loosely model the 155G resonance in 
®^Rb. Explicitly allowing A23 7^ A12 gives qualitatively 
similar results in the three-body calculations described 
below. Allowing A13 ^ provides a direct coupling be- 
tween 1 11) and 1 22) states and, as we discuss later, this 
coupling element is expected to limit the lifetime of Efi- 



mov states near E = e. We note that the methods of [15 1 
permit the application of these ideas over a wider energy 
range because they include the long range van der Waals 
physics in the QDT treatment explicitly. 

In the three-body system, we transform the mass- 
scaled Jacobi coordinates: y^^ = {fj — rk)/d and y^*^ = 
- n) (with d = 21/2/3^/") into hyperspherical 
coordinates by defining the hyperangle a as tana^'-* = 
lyj"'^ l/lyj'''' I, the hyperradius R as = 2/i + 2/2 (which 
is invariant under permutations) and the spherical polar 
angles LOi = {Oi,(j)i\ point in the direction of y^. The 
angular coordinates {a,uji,oj2} are collectively denoted 
f2. The pair-wise interactions are written in the basis of 
two-body internal states, with the third spectator par- 
ticle (near any chosen two-particle coalescence point) in 
either |1) or |2). Hence, if particle 1 is the spectator parti- 
cle, the basis of internal states is {|S)} = {|111), (|112) + 
|121))/\/2, |211),|122),(|212) + |221))/^/2, |222)} with 
energies {Ey,} = {0, e, e, 2e, 2e, 3e}. 



In the multichannel generalization [llj of the adiabatic 
hyperspherical method [17| . we seek solutions to: 
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v{R, n) $(i?; n) = [/(i?)$(i?; n) 
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Here, E_^f^ is a diagonal matrix [i?t/i]ss' = <5ss'£^s, and 
V_{R, 51) is the sum of matrices for each pair- wise inter- 
action expressed in the basis {|E)}. The three-body re- 



duced mass is 11 = m/\/3. The adiabatic potential is 
written in terms of the eigenvalue v = v\ s& U{R) = 
i^i(i^i+4)/2^i?2, where ly^iv^ + A) = 2nR^U{R) - E^). 
Components of the free-space (diagonal) hyperangular 
Green's function G(fJ, il') satisfy, 

(A^ - {ly^ + 4)) Ges (fl, fl') = <5(r!, n') (3) 



Defining a') = N^^i,ijlJ^^_^{a<)fl'^>^^^_{ay), the 

most useful solution is written in the Sturm-Liouville 
form [l^, 

ll .I2 :'"^2 

Yl,^, iuJl)Yi*„,^ {uj[)Yl^,n2 {uJ2)Yil,n2 (^2)] (4) 

where (/^) is a solution to the homogeneous version 
of Eq. ([3]) regular at a = (a = 7r/2), and the normaliza- 
tion Ny^i^i^ is fixed by the Wronskian of and [l^. 
The hyperangular Lippmann-Schwinger (L-S) equation, 



*E(f7) = -~2nR 
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is then solved by evaluating the integral over V^^, in the 
coordinate system where \'lfi'^\ o: \fi — rj\ cx i?sinQ;'^'^^. 
The zero-range s-wave interaction immediately gives li = 
0, and I2 = L. For this Letter, we confine our study 
to states with total orbital angular momentum L ~ 0. 
For equal-mass particles, in the limit a^'^^ ^ we have 
a*-') = a'-^-' = 7r/3; the L-S equation reduces to a matrix 
equation of the form |ll| : 
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where, for bosons. 



A'^^,Xs cot(AE7r/2) 

A(i) -4sin(Ai:7r/6) 
, SE' y3sin(As7r/2) 



= 
(6) 

(7) 

The vector associated with the n'^ eigenstate, Cn"^ in 
Eq. ([6]), is defined by the boundary and normalization 
conditions: 



i = 1 
i = 2,3 



{k} 

E,n 



lim 



g(rfe$E.n) 

drk 



and ^($E,„|$E,„) = 1. (8) 



We have rewritten the eigenvalue a.s = As — 2 
purely for convenience, and P+ and P_ perform cyclic 
and anticyclic permutations respectively upon the basis 
{|E)}. Zeros of the determinant of the matrix in Eq. ([6]) 
yield the eigenvalues A|_„ = 2fj.R'^{Un{R) - Et) + 1/4. 
The potentials Un{R) appear in radial equations of the 
form (ignoring nonadiabatic effects): —F"{R)/2^R^ + 
Un{R)Fn{R) = EF^{R). 
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FIG. 2: (color online) The adiabatic eigenvalues Af (A2) 
with respect to the first (second) threshold are shown for 
a near-resonant choice of e. From left to right, the verti- 
cal dotted magenta lines mark A33, re = 132.284^133, and 
dopen ~ —1.445 X 10*^33. The horizontal dashed red line 
corresponds to the universal eigenvalue for three equal-mass 
resonantly interacting bosons A? = -1.012514 [1, while the 
green dash-dot-dot line indicates the universal eigenvalue for 
only one resonant interaction A^ = 1, and the horizontal dash- 
dot blue line marks the lowest large R free-space eigenvalue 
A^ = 4. 



In Fig. [21 we show eigenvalues with respect to the 
E = and E = e thresholds for e chosen so that a 
quasi-bound two-body state is nearly degenerate with 
the E = Q threshold, giving an identical boson (open- 
open = -1.445 X IOM33. 



channel) scattering length of a, 



Note that we now drop the index n from the label of A 
since the curves shown in Fig. [2] are actually comprised of 
many distinct eigenvalues connected by a series of close 
avoided crossings. This figure illustrates how the mul- 
tichannel problem naturally introduces different length 
scales and symmetry constraints to the three-body prob- 
lem. The eigenvalues may be understood by considering 
simultaneously the internal states of the atoms and the 
important natural length scales: the open-channel 

effective range rg, and the open-channel scattering length 
Oopen • When R <^ A33 , all length scales in our model ap- 
pear long, so with respect to the lowest threshold, we 
have Xf — > —1.012515 appropriate for three resonantly 
interacting bosons. The region A33 <C i? <C is a transi- 
tional regime indicating the importance of the length- 
scale. In the region re <C i? ^ |aopen|, the potential near 
E = behaves again as that of three resonantly interact- 
ing bosons, while near the second threshold E = e, only 
one of the interactions appears resonant and we obtain 



both the E = e and E = 2e thresholds under these con- 
ditions. In Fig. [Sja) we show the lowest 300 adiabatic 
potential curves up to 3.5e. Note the collection of three- 
body curves converging to each three-body threshold at 
E = E^. Note also the quasi-bound two-body thresh- 
old near 0.6e. A three-body collision near this energy 
will couple strongly to the two-body Feshbach resonance 
through this series of avoided crossings. In Fig. Eljc), 
the eigenvalue A2 = A^ — 2fj,R^e is shown to converge 
to the universal value A| = —0.171145 predicted by the 
purely imaginary root of the corresponding single chan- 
nel equation in the limit |a| — > 00 |20| appropriate for two 
identical bosons interacting resonantly with a third dis- 
tinguishable particle of equal mass, consistent with the 
fact that the E ^ e threshold corresponds to internal 
states of the form |E = 2) = (|112) - |121))/^/2 and 
jS = 3) = |211). The spacing of Efimov states obeys 
the formula En = Eq exp {—2Trn / sq) , where A^ — > — s§ 
in the universal regime, and Eq is the energy of the low- 
est Efimov state. If we simply consider the value of the 
Efimov potential at i? ^ 100^433 to be a reasonable esti- 
mate for the energy of the first Efimov state, we obtain 
Eo/h ^ {e/h) — (lOkHz). This places the second state at 
approximately E2/h ^ (e/^) — (ImHz), giving a separa- 
tion of roughly 10 kHz, within the resolution of modern 
radio-frequency experiments [l^ . The spacing between 
states could be made more favorable by using heteronu- 
clear mixtures of atoms, cap italizing on mass-ratios dif- 
ferent from unity. [ll,0,[21|. 

In order for the Efimov states to persist at long- 
range, it is also necessary to show that the positive- 
definite non- adiabatic diagonal coupling —Q(R) = 
— ($(i?)|^3-|$(i?)) corresponding to the Efimov poten- 
tial falls off faster than R^^. This is a non-trivial task 
since the potential is comprised of all potential curves 
approaching the lowest three-body threshold. We de- 
rive an expression for —Q{R) [l8| via a multichannel 
generalization of the method used in [T^]. Our calcula- 
tions indicate that the vector C^'^-' (normalized according 
to Eq. ([8])) for the eigenvalue corresponding to the Efi- 
mov potential has a dominant component: C^}^^, which 
grows linearly with R. Extrapolating this behavior for 
C^i) and using the fact that Xj -0.171145 as R'^, 
the expression for —Q{R) reduces to the simple form: 



the appropriate value A^ = 1 [H. For \aopen\ < R, the -Q{R) = S^^^V^ tt'^^^^, where K4 is the slope of the 



eigenvalues approach the free-space value Af = Aj = 4. 
Note that for values of R where Xf < (Aj < 0) and 
is approximately constant, the potential with respect to 
the E — (E = e) threshold is supercritical and supports 
Efimov states. 

As mentioned before, by considering e to vary with 
the magnetic field, it is possible to make a two-body 
state degenerate with the |cr = 2) threshold. This oc- 
curs when e — ^ l/{mA^^) as illustrated in Fig. [TJ Pro- 
vided that |j4i3| <^ an Efimov potential appears at 



linearly increasing component cj^^\ Since —Q{R) falls 
off faster than R^^, we expect these states to persist at 
long range. 

The lifetime of Efimov trimers near E = e is limited 
by inelastic processes at short and long range. Using a 
WKB estimate, if all inelastic transitions occur at short 
range with probability Psr (a reasonable assumption in 
the limit A13 0), then the width of successive reso- 
nances will scale geometrically with the n}^ energy level 
as TsR « Psfl^oe-2""/^Vso. 
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FIG. 3: In (a), we show approximately the lowest 300 potential curves up to 3.5e, while (b) shows an enlarged view of the 
region near the E = t threshold, and the attractive diabatic Efimov potential. Note also the series of avoided crossings in 
(a) near E = 0.6e indicating the presence of a two-body quasi-bound state. In (c) we show the eigenvalue near E = e converging 
to the universal value for two identical bosons and one distinguishable particle A'^ —0.171145 indicated by the dashed red 
line. 



For finite Ai^ inelastic transitions can also occur at 
long range due to two-body inelastic transitions, and the 
width of Efimov resonances is dominated by the size of 
Ai3. Letting A13 ^ serves to broaden the avoided 
crossing comprising the universal Efimov potential in 
Fig.EJb), dramatically increasing the decay rate to three- 
body continuum states associated with the |S = 1) state, 
and making the Efimov states shorter lived and more dif- 
ficult to observe. This is understood by noting that when 
the coupling between |cr = 3) and \a = 1) is significant 
(i.e. Ayi is sizable), direct inelastic transitions between 
|(T = 3) and \a = 1) states dominate over inelastic tran- 
sitions via the intermediate \a ~ 2) state. 

To summarize, we have identified a new class of Efimov 
states embedded in the three-body continuum. The hy- 
perradial potentials supporting these states have univer- 
sal properties consistent with the symmetry constraints 
implied by the internal degrees of freedom of the three- 
atom system, and can further be understood in terms of 
relevant two-body length scales. Further, we have identi- 
fied the coupling constant A13 responsible for limiting the 
lifetime of these Efimov states at long range, and have 
shown that the repulsive diagonal nonadiabatic correc- 
tion falls off as i?^'^, consistent with the single-channel re- 
sult. We stress that since the potential supporting these 
states appears when a quasi-bound two-body state is de- 
generate with an excited two-body threshold, the open- 
channel scattering length in general will not be resonant, 
and the gas is expected to be stable with respect to the 

recombination scaling law. We postulate that it may 
be possible to observe these states by spectroscopic tech- 
niques, perhaps with sufficient accuracy to measure two 
Efimov resonances for the first time. 
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